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A  HEURISTIC  METHOD  OF  OPTIMAL  GENERALIZED 
HYPERCUBE  ENCODING  FOR  PICTORIAL  DATABASES 


INTRODUCTION 

In  similarity,  retrieval  from  a  pictorial  database  [1] ,  pattern  recognition  [2,4] ,  and 
clustering  analysis,  it  is  often  desired  to  find  the  set  of  database  records  (or  set  of  patterns, 
n-dimensional  feature  vectors,  etc.)  that  are  most  similar  to  a  test  record  (or  test  pattern, 
test  feature  vector,  etc.).  In  the  case  of  large  databases  (or  patterns,  clusters),  it  is  very 
useful  to  encode  the  original  database  into  certain  convenient  format  in  order  to  facilitate 
similarity  retrieval  and  updating. 

In  this  report,  we  shall  consider  one  such  technique  called  “Generalized  Hypercube” 
encoding  [3]  and  describe  a  heuristic  method  of  generating  minimum  number  of  GH 
encoded  tuples. 


GENERALIZED  HYPERCUBE  ENCODING 

Given  a  set  of  points  S  =  {  p* ,  p2,  • .  pk  }  >n  n-dimensional  space,  where  ps  = 

(xx ,  x2, . .  .,  xn)  denotes  a  point  in  S,  we  can  have  a  family  of  n-dimensional  GH  codes. 
For  each  m,  1  <L  m  n  +  1,  the  GHm  codes  are: 

(Xj , . . xm_j ;  am ,  .  .  .,  a,, ;  bm, . . bn). 


where  (Xj ,  x2, .  .  .,  xm  l ,  z  , . .  .,  zn)  is  in  S  for  some  coordinates,  zm , .  .  .,  zn;  and 

a. j  =  min  {y^  for  some  coordinates,  zk, (xx , .  .  .,  x^ ,  zm, ....  yj, ....  z„)  in  S  }  , 

b.  =  max  {  y}:  for  some  coordinates,  zk,  (xv  . . .,  x^,  zm, . . .,  . . zj  in  S}  . 

When  m  =  1,  we  are  simply  using  the  smallest  n-dimensional  hypercube  containing  S 
as  the  GH  code.  When  m  =  n+1,  the  original  point  set  S  is  used  as  the  GH  code.  Other  in- 
between  values  of  m  give  GH  codes  of  various  levels  of  details. 


For  example,  if  point  set  S  =  { (1,1, 3),(1, 1,5), (1,2,1), (1,2,4), (2,3,6) }  ,  GHj  is 
{ (1,1,1;2,3,6)}  .  GH2  is  { (1;1,1;2,5),(2;3,6;3,6) }  .  GH3  is  { (1,1;3;5),(1,2;1;4),(2,3;6;6)  }  . 
GH4  is  S  itself.  It  should  be  noted  that  other  GH  codes  can  be  generated,  if  we  permute  the 
coordinates.  The  previous  example  we  have  {(3;1,1;1,1),(5;1,1;,1,1),  (1;1,2;1,2), 

(6;2,3;2,3) }  . 


Therefore,  it  is  important  to  select  carefully  the  coordinates  on  which  the  GH  encoding 
technique  can  be  applied,  i.e.  to  generate  a  minimum  number  of  GH  encoded  tuples  the 
problem  reduces  to  select  a  special  "handle”  set.  It  is  easy  to  see  that  once  a  handle  vector 
has  been  selected,  the  GH  encoding  is  unique. 
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In  this  report  we  shall  describe  the  following  heuristic  methods  for  problems: 

•  Select  a  handle  set  {ij,  i2, . . }  from  { 1,  2, . .  n  }  so  that  the  least  number 
of  GH  encoded  tuples  are  generated. 

•  Suppose  we  define  the  density  of  a  hypercube  H  as  the  number  of  data  points  in  H, 
divided  by  total  number  of  lattice  integers.).  For  example  the  hypercube  (1,1,1  ;2,3,6) 
previously  described  has  a  density  of  5/2 *3*6  =  0.1388.  Intuitively,  we  would  like  to 
generate  hypercubes  having  high  densities.  Therefore  Problem  2  can  be  stated  as  follows: 

Given  a  threshold  density  t,  0  t  <.  1 .0,  find  optimal  GH  encoding  such  that  each 
hypercube  has  density  no  less  than  t. 


PROBLEM  1 

Given  m,  the  method  of  choosing  m-1  handles  which  will  generate  optimal  GHm  codes 
for  a  point  set  S  is  based  on  the  following  ideas: 

•  For  each  handle  (column)  i,  (i  =  1, 2, . . .,  n)  consider  vector  Vj  =  (xli,  x2j, .  .  .,  xkj) 
where  k  is  the  number  of  points  in  the  set  S.  Among  the  elements  x;i ,  x-2, .  . .,  xik  those 
having  the  same  value  are  grouped  together  and  their  associated  points  (in  S)  will  consist 

of  disjoint  groups,  say,  Gjj,  G2i, . . .,  G8i.  Define  the  count  (or  cardinality)  of  G.(  =  C.(i). 

We  next  compute  a  measure  of  the  priority  of  the  ith  column  p..  ‘ 

fi 

•  Pi  -  2  rcj(i)-112- 

j=i 

Intuitively,  we  ignore  those  groups  having  only  one  element  and  put  heavier  weight  to 
those  groups  having  more  than  one  element.  Obviously,  there  are  many  other  ways  to  define 
a  priority  measurement  for  each  handle. 

•  Based  on  the  measure  P.,  (i  =  1,  2, . . .,  n)  choose  the  best  m-1  handles,  i.e.,  i’s  which 
have  larger  values  of  Pj. 

Following  is  the  algorithm  to  implement  the  above  ideas: 


ALGORITHM 
Step  1. 

To  conpute  P(,  consider  vector  V.  =  (Xj.,  x2i, . . .,  xkj), 
a.  Make  a  stack  of  all  k  elements  of  V. . 


,a3rtL^-»  - - - 
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b.  Set  a  test  element,  TELMNT  =  top  element  of  stack  and  set  COUNT  =  1 . 

c.  Replace  the  top  element  of  the  stack  by  the  last  element  and,  set  k  =  k-1 ;  if  k  <  1, 
go  to  substep  e. 

d.  Restack  the  k  elements  to  check  that  the  top  element  =  TELMNT. 

•  If  it  does,  COUNT  =  COUNT  +  1;  then  go  to  substep  c. 

•  If  it  does  not,  then 

{  If  COUNT  >  1,  then  P;  =  ?i  +  COUNT  *  COUNT  and  COUNT  =  1.  TELMNT  = 
top  element  of  the  stack;  go  to  substep  c  } . 

e.  Include  the  last  group’s  population  count  in  Pj 

If  COUNT  >  1,  then  R  =  P;  +  COUNT  *  COUNT. 


Step  2. 

Using  STEP  1,  compute  all  P.’s  for  i  =  1,  2,  . .  .,  n.  Store  handle  numbers  in  an  array, 
such  as  [c(i)  =  i,  i  =  1, . .  .,  n] . 

Step  3. 

Sort  array  P  in  descending  order,  and  while  moving  the  elements  of  P  during  sort,  move 
elements  of  c  also  in  the  same  fashion. 


For  generating  GH  codes,  choose  handles  in  the  first  (m  -  1)  positions  of  array  c, 
which  correspond  to  the  best  (m  -  1)  values  of  Pj,  i  =  1,  2, . .  .,  n. 

Once  the  set  of  handles  has  been  selected,  the  GH  encoding  is  unique. 


For  example,  let  point  set  S  =  { (3,1,1),  (5,1,1),  (1,1,2),  (4,1,2),  (6,2,3) }  .  Using  the 
previous  algorithm,  the  measure  of  priority,  Pj  for  i  =  1,2,3  are:  Pj  =  0,  P2  =  16,  P3  =  8. 
Based  on  the  best  values  of  P. ,  we  choose  handle  set  (x2  ,x3  ,x. )  and  with  respect  to  this  set 
of  handles  (GH.)  is  { (1,1,1;2,3,6) }  .  GH2  is  { (1;1,1;2,5),  (2;3,6;6) } .  GH,  is  { (1,1;3;5), 
(1,2;1;4),(2,3;6;6)}  .GH4  is  { (1,1,3),  (1,1,5),  (1,2,1),  (1,2,4),  (2,3,6) }  . 


EXPERIMENTAL  RESULTS 

The  previous  algorithm  was  implemented  on  an  IBM/370  in  FORTRAN;  a  listing  of  the 
program  is  given  in  Appendix  A.  The  data  for  the  program  were  random  numbers  with 
normal  distribution  and  standard  deviation,  a  number  for  each  column  vector  V.  was  also  a 
random  number  with  uniform  distribution. 


I 
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The  selected  set  of  handles  by  the  previous  method  was  used  to  generate  GHm  encoded 
tuples,  and  their  number  was  counted  for  a  different  number  of  point  set  and  dimensions 
of  space.  To  find  the  optimal  GHm  encoding,  great  effort  was  taken,  and  the  GH  codes 
with  minimum  number  of  GHm  tuples  was  taken  as  the  optimal  GHm  encoding.  (Bee  Tables 
1,  2,  and  3.) 

Few  other  heuristic  approaches  were  tested  to  select  the  proper  set  of  handles  to 
generate  optimal  GH  encoding,  but  the  previous  method  has  shown  the  best  results  and  will 
be  used  in  the  next  optimization  problem. 


PROBLEM  2 

From  the  experimental  results  we  know  that  the  smaller  the  value  of  m  (length  of 
handles)  the  smaller  the  number  of  GHm  encoded  tuples  for  a  point  set  S,  but  at  the  same 
time  such  hypercubes  tend  to  be  sparse.  Therefore,  to  optimize  GH  codes  for  a  given 


Table  1  —  The  number  of  optimal  GHm  encoded  tuples 


No.  of 
Points 
(M) 

gh2 

gh3 

H 

0 

Error 

(%) 

H 

O 

Error 

(%) 

10 

8 

8 

0.0 

10 

10 

0.0 

20 

15 

15 

0.0 

20 

20 

0.0 

30 

22 

22 

0.0 

30 

30 

0.0 

40 

12 

12 

0.0 

33 

33 

0.0 

50 

34 

34 

0.0 

50 

50 

0.0 

60 

16 

16 

0.0 

56 

55 

1.02 

70 

24 

24 

0.0 

68 

68 

0.0 

80 

30 

30 

0.0 

75 

75 

0.0 

90 

27 

27 

0.0 

88 

87 

1.16 

100 

20 

20 

0.0 

86 

86 

0.0 

Dimension  N  ■  3 
H  ■  heuristic 
O  •  optimum 
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Table  2  —  The  number  of  optimal  GHm  encoded  tuples 


No.  of 

gh2 

a 

jr 

CO 

gh4 

Points 

(M) 

H 

O 

Error 

(%) 

H 

O 

Error 

(%) 

H 

O 

^  ■ 
El 

10 

8 

8 

0.0 

10 

10 

0.0 

10 

10 

0.0 

20 

16 

16 

0.0 

20 

20 

0.0 

20 

20 

0.0 

30 

10 

10 

0.0 

29 

28 

3.57 

30 

30 

0.0 

40 

21 

21 

0.0 

40 

39 

2.56 

40 

40 

0.0 

50 

16 

16 

0.0 

44 

44 

0.0 

50 

50 

0.0 

60 

25 

25 

0.0 

59 

59 

0.0 

60 

60 

0.0 

70 

28 

28 

0.0 

68 

68 

0.0 

70 

70 

0.0 

80 

21 

21 

0.0 

74 

71 

4.23 

80 

80 

0.0 

90 

18 

18 

0.0 

75 

75 

0.0 

90 

90 

0.0 

100 

22 

22 

0.0 

93 

93 

0.0 

100 

100 

0.0 

Dimension  N  *  4 
H  «  heuristic 
O  »  optimum 
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Table  3  —  The  number  of  optimal  GHm  encoded  tuples 


No.  of 

gh2 

gh3 

gh4 

gh5 

Points 

(M) 

H 

O 

Error 

(%) 

H 

O 

Error 

(%) 

H 

O 

Error 

(%) 

H 

O 

Error 

(%) 

10 

8 

8 

0.0 

10 

10 

0.0 

10 

10 

0.0 

10 

10 

0.0 

20 

13 

13 

0.0 

20 

19 

5.26 

20 

20 

0.0 

20 

20 

0.0 

30 

10 

10 

0.0 

29 

29 

0.0 

30 

30 

0.0 

30 

30 

0.0 

40 

16 

16 

0.0 

38 

38 

0.0 

40 

40 

0.0 

40 

40 

0.0 

50 

21 

21 

0.0 

49 

49 

0.0 

50 

50 

0.0 

50 

50 

0.0 

60 

18 

18 

0.0 

55 

55 

0.0 

59 

59 

0.0 

60 

60 

0.0 

70 

18 

18 

0.0 

64 

63 

1.58 

70 

70 

0.0 

70 

70 

0.0 

80 

21 

21 

0.0 

75 

75 

0.0 

80 

79 

1.26 

80 

80 

0.0 

90 

26 

26 

0.0 

84 

84 

0.0 

90 

90 

0.0 

90 

90 

0.0 

100 

22 

22 

0.0 

96 

92 

4.35 

100 

100 

0.0 

100 

100 

0.0 

Dimension  N  *=  5 
H  =  heuristic 
O  =  optimum 


threshold  density  t,  find  a  set  of  handles  with  minimum  length  m,  so  that  the  resultant 
hypercubes  have  densities  no  less  than  theshold  density,  t. 

Two  special  cases  can  be  solved  immediately.  If  t  =  0,  then  we  can  select  m  =  1,  so  that 
a  single  n-dimensional  hypercube  is  generated,  and  GHj  becomes  the  optimum  GH  en¬ 
coding.  If  t  =  1.0,  then  all  the  hypercubes  must  have  unit  density,  and  we  can  select 
m  =  n  +  1  so  that  the  original  point  set  S  becomes  the  optimum  GH  encoding.  For  other 
values  of  t,  m  should  fall  between  1  and  n  +  1. 
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ALGORITHM 

The  algorithm  optimize  GH  encoding  for  a  given  threshold  density. 

Step  1. 

Start  with  m  =  1; 

generate  GHX  codes,  which  is  one  hypercube.  If  the  density  of  this  hypercube  is 
>  =  threshold  density,  then  this  is  the  optimum  GH  encoding,  and  go  to  Step  7, 

otherwise,  set  m  =  m  +  1. 

Step  2. 

By  using  the  previous  heuristic  method  discussed  in  Problem  1,  choose  a  set  of  handles 
which  will  generate  the  optimal  GHm  encoding  for  a  given  m. 

Set  a  flag  =  0;  and 

NUM  =  0,  where  NUM  will  be  the  number  of  GHm  encoded  tuples  with  respect  to  the 
chosen  handle  vector  of  length  m. 

Step  3. 

a.  Generate  the  next  hypercube  with  respect  to  the  above  chosen  handle  vector  and 
compute  its  density. 

b.  If  density  of  this  hypercube  =  >  threshold  density,  and 
if  the  number  of  points  in  this  hypercube  >  1,  then 

{  include  this  hypercube  in  optimal  GH  encoding. 

NUM  -  NUM  +  1;  and 

delete  the  points  included  in  this  hypercube  from  the  original  point  set.  } 

c.  If  the  density  of  this  hypercube  <  threshold  density,  then  set  the  flat  =  1 
Step  4. 

Repeat  Step  3  until  all  the  points  in  point  set  have  been  considered  once. 

Step  5. 

If  the  number  of  points  in  remaining  point  set  =  0  then  go  to  Step  7. 

If  (number  of  points  in  point  set  =  1  or 
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NUM  <  =  0  and  flag  =  0),  then  go  to  Step  6. 

If  (NUM  <  =  0  and  flag  =  1),  then  m  =  m  +  1. 

Go  to  Step  2. 

Step  6. 

Generate  GHm  encoded  tuples  with  respect  to  the  last  handle  vector  of  length  m. 
Each  GHm  encoded  tuple  is  included  in  the  optimal  GH  encoding. 

Step  7. 

Stop. 

As  an  illustration,  if  point  set  S  =  {  (3, 1,1), (5, 1,1), (1,1, 2), <4, 1,2), (6, 2, 3)  }  and 
threshold  density  =  50%  (0.5),  then  with  the  previous  algorithm,  optimal  encoding  is 
attended  as  { ( 1  ;1,2;2,5),  (2;3,6;3,6) }  and  the  density  of  each  hypercube  is  >  =  50%. 


EXPERIMENTAL  RESULTS 

The  previous  algorithm  also  was  implemented  on  an  IBM/370  in  FORTRAN.  A  listing 
of  the  program  is  given  in  Appendix  B.  The  test  data  were  again  random  numbers. 

The  results  are  in  Tables  4,  5,  and  6. 


DISCUSSION 

In  this  section,  some  interesting  applications  of  the  generalized  GH  encoding  technique 
are  discussed. 

Reference  1  suggests  that  a  relational  file  R  can  be  characterized  by  a  logical  ex¬ 
pression  E,  so  that  every  n-tuple  (xj ,x2, . . .,  xn)  in  R  satisfies  this  logical  expression. 

The  generalized  GH  encoding  technique  can  be  used  to  find  a  suitable  logical  ex¬ 
pression  characterizing  a  relational  file. 

For  example,  if  S  is  as  given  in  the  Generalized  Hypercube  Encoding  section  of  this 
report  and  m  =  1,  then  GHj^  =  { (1,1,1;2,3,6) }  and  the  corresponding  logical  expression  is 
(1  <  =  Xj)  &  (1  <  =  x2)  &  (1  <  =  x3)  &  (Xj  <  =  2)  &  (x2  <  =  3)  &  (x3  <  =  6); 

when  m  =  2,  then  GH2  =  { (1;1,1;2,5),(2;3,6;3,6) }  and  the  corresponding  logical 
expression  is  (Xj  =  1)  &  (1  <  =  x2)  &  (1  <  =  x3)  &  (x3  <  =  5))  V  ((xa  =  2)  &  (x2  =  3) 

&  (x3  =  6)),  thus,  each  GH  encoding  corresponds  to  a  logical  characteristic  expression  for 
the  relational  file. 
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Table  4  —  The  number  of  optimal  GH  encoded  tuples 


(M) 

t=10 

t=20 

t=30 

t=40 

t=50 

t=60 

t=70 

t=80 

t=90 

10 

6 

7 

7 

7 

9 

9 

10 

10 

10 

20 

11 

14 

18 

18 

18 

18 

18 

18 

18 

30 

9 

23 

24 

24 

24 

26 

28 

28 

28 

40 

18 

36 

36 

36 

36 

37 

38 

38 

38 

50 

23 

36 

39 

40 

42 

48 

50 

50 

50 

60 

12 

43 

43 

48 

53 

55 

58 

58 

58 

70 

mm 

40 

46 

46 

49 

52 

64 

64 

64 

80 

H 

52 

53 

54 

57 

57 

60 

60 

60 

90 

45 

50 

52 

56 

77 

80 

84 

84 

100 

■ 

40 

52 

58 

60 

81 

82 

91 

91 

M  <■  number  of  points  in  point  set 
t  *  percent  threshold  density 
Dimension  N  “  3 


9 


YANG,  CHANG  AND  SINGH 


Table  5  —  The  number  of  optimal  GH  encoded  tuples 


(M) 

t=10 

t=20 

t=30 

t=40 

t=50 

t=60 

t=70 

10 

7 

7 

9 

9 

9 

9 

10 

20 

14 

17 

17 

18 

18 

18 

18 

30 

25 

28 

28 

29 

29 

29 

29 

40 

32 

36 

39 

39 

39 

39 

39 

50 

38 

45 

46 

49 

50 

50 

50 

60 

46 

52 

56 

57 

57 

57 

59 

70 

42 

57 

60 

68 

68 

68 

70 

80 

45 

72 

77 

78 

79 

79 

79 

90 

58 

72 

78 

79 

82 

87 

88 

100 

66 

78 

81 

90 

91 

92 

_ 

92 

M  =  number  of  points  in  point  set 
t  -  percent  threshold  density 
Dimension  N  *  4 


i 
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Table  6  —  The  number  of  optimal  GH  encoded  tuples 


(M) 

t=10 

t=20 

t=30 

t=40 

t=50 

t=60 

t=70 

t=80 

t=90 

10 

10 

10 

10 

10 

10 

10 

10 

10 

10 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

30 

28 

29 

30 

30 

30 

30 

30 

30 

30 

40 

40 

40 

40 

40 

40 

40 

40 

40 

40 

50 

47 

50 

50 

50 

50 

50 

50 

50 

50 

60 

57 

58 

58 

58 

58 

59 

59 

59 

59 

70 

68 

68 

69 

69 

69 

70 

70 

70 

70 

80 

76 

77 

77 

77 

77 

79 

79 

79 

79 

90 

79 

86 

87 

87 

88 

89 

90 

90 

90 

100 

95 

100 

100 

100 

100 

100 

100 

100 

100 

M  “  number  of  point*  in  point 
t  “  percent  threibold  density 
Dimension  N  -  6 


In  addition  to  applications  in  database  characterization,  the  previously  described 
technique  also  can  be  applied  in  clustering  analysis  where  clusters  are  described  by  hy¬ 
percubes  with  handles.  Thus  optimization  technique  will  generate  optimal  encoding  for  the 
clusters. 

A  third  application  can  be  the  description  of  n-dimensional  pictures,  which  will  have 
a  practical  significance  in  computer  graphics  and  image  processing  applications. 
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APPENDIX  A 


1 .  //JUM<  J33 

2.  T=( 8.0  I .L=2.R=2S6 

3.  //  TXEC  FORTGCLG 

4.  //^OPT  .SySIN  DO  * 

5.  CCCCCCCCCC 


6. 

c 

7  . 

c 

MAIN  PR1GRAM  TO  GPNFR4TE  OP  T I M4L  G«M  FNCODING 

FC» 

4  POINTSET  BY 

P. 

c 

A  HEJPISTIC  METHOD  t  TO  CHECK  T HF  RESULTS.  NUMBER 

OF  THE  SF 

GHM 

9. 

c 

FNCDjED  *U3LF4  IS  COMPARED  * ITH  the  MINIMUM  (REAL 

OPT IMUN  » 

NO  • 

10  . 

c 

OF  GrIM  ENCODEO  TUPLES  FOUND  IN  EXHUASTIVF  APPROACH. 

1 1  . 

c 

1  2  . 

c 

cal.s; 

1  3. 

c 

(11  GENPTS:  TO  GENERATE  TEST  DATA  I.E.  A  POINT 

SET. 

1  4  . 

c 

(2)  hfuris:  to  generate  optimal  ghm  encoded 

TUPLES  BY  A 

15. 

c 

HEURISTIC  METHOD. 

16, 

c 

( 3 1  EXNAUS:  TO  GENERATE  OPTIMAL  GHM  ENCODING 

.  BY 

EXHAUSTIVE 

17. 

c 

APPROACH • 

IP. 

c 

19. 

ccecccrcce 

20  . 

integer  ptseti i i o.6» .pi 6i .ncdde (5.31 

21  . 

TEL=1 

22. 

N  =  4 

23. 

M  =  50 

2*. 

CALL  GENPTSIPTSET.M.NI 

’5. 

IM ICL  .EO.  11  GO  TO  440 

26. 

PRINT  1 

27. 

I 

0  PM  A  T I  • 1  1  « 1 0  X  .  •  INIT! AIL  POINT  SET  1S:»> 

75. 

D3  100  I  *  l .M 

29. 

PRINT  2  . (PTSFTI I , Jl . J=l ,NI 

30. 

1  00 

CONTJ  NUE 

31  . 

2 

FORMAT! '0* . I0X  ,• (  •  .513. •  1 • I 

32. 

PRINT  3 

3. 

3 

=0RMAT(////, 1001 •♦•II 

1  . 

PRINT  4 

35. 

4 

■  3  R  m  A  T !// . 75X  « • HEURISTIC  APPROACH*  1 

36. 

PRINT  5 

37. 

5 

fqpm4T(  J»5X  *  • - - •  1 

35. 

PRINT  6 

39. 

6 

FORMAT!//,  1001  •♦•  1  I 

40. 

440 

CALL  HEUPJ  S(PT  SET .M.N. NCODE ,  IFL 1 

41  • 

!=t I=L  .EQ.  II  GO  TO  442 

4?. 

PRINT  U 

43. 

1  1 

FORMAT {///, 100 (  •  4»  n 

44  • 

PRINT  12 

45. 

1  2 

FORMAT {// %? IX.  *  E  XH4US 7  f  C  APPROACH’  1 

4f>. 

PRINT  13 

47. 

1  3 

FORMAT!  2  IX  ,  • - - — - •  1 

4  R  . 

PRINT  14 

49  • 

1  4 

FORMAT  I //,  100!  •♦  •  1  1 

50  - 

44? 

CALL  F  XH AUSIPTSET .M.N. NCODE. IFL 1 

51  . 

Nl  =  N-l 

52  . 

PRINT  20. M 

53. 

20 

F(jPMAT  I//.I0X,  'NO.  OF  POINTS  IS  *  ••  141 

54. 

PRINT  22. N 

55. 

22 

format ! i ox . »no .  of  oimensins  *  *.  !3i 

56. 

PRINT  15 

57  . 

1  5 

FORMAT  I / / . 20  X •  'HEURISTIC*  . 1 0 X . • OPT  I  MUM •  ,5X , 

*  WORST '  *  BX , *  KERR 1 

PR. 

DJ  200  I  =  I.NI 

5°. 

ERR  =  (NCODE!  I.  II  -  NC ODF (  I . 2  I  > ♦ 1 00. /NCODE ! 1 . 2 1 

60  * 

2  00 

PRINT  16.  ! •  (NCOOE(!.JI.Jb  1.31. FRR 

13 
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61  . 

16 

=  3PM*T(I0'  . I0X. • GH*  . 12.3(1  OX. 141. F16. 31 

62  . 

IF( IFL  .50.  1  1  POT  NT  66 

63. 

66 

=ORMAT1//,IOOI  •  .*  I 1 

6*  . 

500 

CONTI NUE 

65. 

600 

CONTI NUE 

66. 

STOP 

67. 

END 

6  0. 

CCCCCCCCC 

69. 

C 

70. 

c 

R4ND3«(Xt:  GENERATES  A  RANDOM  NUMBER  X  •  0<=  X 

<=1.0  WITH 

71  . 

c 

UNIFORM  DISTRIBUTION* 

72. 

c 

7?  . 

cccccccccc 

74. 

suanouTiNE  randomixi 

75. 

REAL  X 

76. 

INTEC.EP  A/19727/.  MULT/2521  1/.  BASE/32760/ 

77. 

A  =  M3DIMULT4 A. BASE t 

70. 

X  =  F.  0  AT 1  A  1/FLOAT ( BASE  1 

79. 

RETURN 

no. 

end 

01  . 

cccccccccc 

82. 

c 

83  . 

c 

genpts(ptset,m,ni to  generate  test  data  set 

IN  APR  AY 

PTSET. 

84  . 

c 

•MICH  contains  m  number  of  points  in  n  dimensional 

SPACE. 

85. 

c 

calls: 

06. 

c 

(II  FA0345  A  LIBRARY  SUBROUTINE  TO  GENERATE 

RANDOM  NUMBERS 

07. 

c 

NORMAL  DISTRIBUTION  WITH  A  GIVEN 

standard 

88. 

c 

DEVIATION. 

09. 

c 

(2i  random:  to  get  anothfr  random  number  to 

BE  USED 

AS  FIX 

90. 

c 

STANDARD  DEVIATION  FOR  ONE  COLUMN. 

6*  . 

c 

• 

cccccccccc 

-  J  • 

SJBROJTINE  GENPTSIPTSET.  M.NI 

94. 

INTEGER  PTSET ( 110.61 

95. 

03  100  J  s  I.  N 

96. 

CALL  R  ANOCM ( U • 

97. 

STD  *  I9.04U  ♦  3 

95. 

93  100  I  «  1,  M 

99. 

CAL.  FA03AISTO.XI 

100. 

PTSFT ( I .J 1  =  INT( X  ♦  .5 1 

101  . 

100 

CONTINUE 

1  02. 

PE  TURN 

1  03. 

END 

104. 

cccccccccc 

1  05. 

c 

1  06. 

c 

pROCES(RTSET.o.m,ni:-  COMPUTES  THE  MEASURE  of 

PRIORITY 

P(  1  1 

1  07. 

c 

FOP  EACH  COLUMN  I  WHERE  I  =  1.....N  FOR  A 

GIVEN  POINT  SET 

1  08. 

c 

IN  ARRAY  PTSET  CONTAINING  M  POINTS  IN  N  DIMENSIONAL  SPACE 

1  09. 

c 

calls: 

I  10. 

c 

(II  HPOOER:  TO  COMPUTE  PHI  FOR  ONE  PARTICULAR  VALUE 

OF  I  . 

111. 

c 

112. 

cccccccccc 

1  13. 

SUaRDUTINE  PROCESI PTSET ,D, M.NI 

114. 

INTEGER  3TSET (  l  10,61 , TE  MPR(  1 101 ,P( 61 

1  15. 

DO  IO0  I  =  l.N 

1  1  6. 

OD  200  J  =  1 .M 

t  17. 

POO 

TEWPBI J  1  =  PTSET( J.I  * 

110. 

CAt_L  HPOOERI  1  .TEMPR.M.Pl 

1  19. 

1  00 

continue 

120. 

RETURN 

14 


121  . 
122. 
123. 
12*. 
I  25. 
I  26. 
127. 
1  26. 
129. 
133. 

131. 

132. 

133. 

134. 
I  35. 

136. 

137. 
136. 
I  39. 
140. 
141  . 

142. 

143. 

144. 

145. 

146. 
1  47. 
146. 

149. 

150. 
51. 

.2. 
i  53  . 

154. 

155. 

156. 

157. 
156. 

159. 

160. 
161  . 
162. 

163. 

164. 
I  65. 
166. 
1  67. 
1  68. 

169. 

170. 
I  7  |  . 
172. 
I  73. 
I  74  . 
1  75. 
I  76. 
177. 
1  76. 
179. 
160. 
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593 

CCCCCCCCCC 

c 

C  MB03ER1 1 . .M .O I S-  COMPUTES  P( I  *  6  09  «  COLUMN  VECTOR  TEMPP 
C  4NO  USES  HE  49  STRUCTURE. 

C  C4t_L  S  : 

C  <11  9£HE4P1  to  MAKE  HE 4 P  06  ELEMENTS  06  4RP4V  TEMPP. 

C 

CCCCCCCCCC 

SUBROUTINE  HPODERI 1 .TEMPP, M. PI 
1 N  TEG: P  TEMPP! 1  l 0  1 .9(6 t 
INTEGER  TEMP. ELMSNT ,PSUM. SUM 
Ml  ■  M 
L  «  M/2  ♦  l 
10  L  «  L-  I 

C4LL  REHEAPIL.MI  .TEMPRI 
!c  <.  .GT.  II  GOTO  10 
elment  =  0 

PSUM  =  0 

SUM  *  0 

20  TEMP  =  TEMPP  f 1  I 

TEMPP  <11  *  Tf  MPP<  Ml  1 

T  6  IE.MENT  .EO.  TEMP1G0  TO  30 
I6ISUM  .GT.  1  I  PSUM  a  PSUM  »  SUM  •  SUM 
ELMENT  a  TEMP 
SUM  =  l 
GO  TO  40 

30  SUM  •  SUM  ♦  1 

40  M|  a  M|  -  I 

I 6 ( M I  .GT.  01  GO  to  50 

16  < SUM  .GT.  It  PSUM  *  PSUM  4  SUM  »SUM 

P<I I  =  PSUM 

RETURN 

50  CALL  PFHE4PI I .Ml .TEMPP I 

GOTO  20 
END 

CCCCCCCCCC 

c 

C  REHEAR!. .  M  I  .TCM9R  •  MAKES  HEAP  06  ELEMENTS  06  TFMPR  STARTING 
C  c ROM  TFMPR<L I  TO  TEMPRIMI). 

C 

CCCCCCCCCC 

SUBROUTINE  REME»P<L.M1 .TEMPRI 
INTEGER  TEMPRI l 101 
INTEGER  FLAG  .  X 

rr  «  l 

JJ  »  2* 1 T 
X  r  TEMPRI  I  I  I 
6_AG  *  1 

16  <JJ  .GT.  Mil  RETURN 
40  |6<JJ  ,GE.  Mil  GOTO  10 

16  I  TEMPT  <  J J  I  .LT.  TEMPRI l 1 1  JJ  aJJ*| 

10  I  6  I  X  . GE .  TEMPP I J  J  t I  GOTO  20 
TEMPRI! II  a  TEMPRI JJ I 
1  I  =  JJ 
JJ  =  2MI 
TEMPRI  11  I  a  X 
GOTO  30 

20  FLAG  =  0 

30  1 6  ( JJ  .LE.  Ml  .AND.  FLAG  .EO.  It  GOTO  40 
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lev . 

RETURN 

IS?. 

END 

1  B3. 

cccccccccc 

IB*. 

r 

IBS. 

c 

SORT! P. 0. N 1 : -  SORTS  P  ARRAY  CONTAINING  N  ELFMFNTS  ANO  DUPING 

1  B6. 

c 

THIS  SORT  ALSO  MOVES  THE  ELCMENTS  OF  ARRAT  0  IN 

THE  SAME 

1ST  . 

c 

FASHION. 

188. 

c 

1  89. 

cccccccccc 

190. 

SU390UT1NE  S0RT1P.0.NI 

191  . 

INTEGER  0(6 I.PI6) 

192. 

INTEGER  FLAG.  TOP.  TEMP.  BOTTOM 

193. 

flag  ■  1 

1  94  * 

BOTTOM  *  N 

195. 

TO*  «  1 

196. 

i  0 

IFIFLAG  .NF.  0  .AND.  TOP  .LT.  BOTTOM  I  GO  TO  20 

197. 

RETURN 

198. 

20 

FLAG  *  0 

109. 

II  *  BOTTOM  -l 

200  . 

00  100  I  *  1.11 

201  . 

if  (P(!l  .r.E.  PI1MII  GO  TO  100 

202. 

*  «  PI  I  1 

20  3. 

Pill  =P( I M  1 

704  . 

P(  1  *1  l  *  * 

205. 

temp  *n( i » 

206. 

01  1  1  =0(1 H  1 

207. 

0( (*1 IxTFMP 

208. 

Flag  =  1 

209. 

1  00 

CONTINUE 

210. 

93TT0M  -  bottom  '1 

'•  1  • 

GO  TO  10 

2. 

END 

-13. 

cccccccccc 

214. 

c 

215. 

c 

soptpt (ptset .o.m.n i : -  t0  sort  points  of  pointset  in 

PTSET  M.R 

216. 

c 

hANOlF  VECTOR  in  ARRAV  0.  HEAP  SORT  IS  USED  FOR 

this. 

217. 

c 

Calls: 

21  8. 

c 

tit  renip:  to  make  heap. 

219. 

c 

(21  C0N9AR:  TO  COMPARE  ORDER  OF  TWO  POINTS. 

270. 

c 

(31  EXCH4N:  TO  EXCHANGE  POSITIONS  OF  TMO  POINTS  in  ptset. 

221  . 

c 

222. 

cccccccccc 

223. 

SUBROUTINE  S0RTPT1PTSET, O.M.N! 

224. 

INTEGER  RTSETI 1 1 0.61 .0(6 1 

225. 

Ml  x  M 

226. 

L  *  M 1  /Z  ♦  1 

227  . 

10 

L  x  L  -  1 

228. 

CALL  REM|P|L,M1 .M.N.O.PTSETI 

229. 

I  c  I  L  .GT . 1  |  GOTO  1 0 

230. 

20 

CALL  F XCHANI 1 .Ml ,N, PTSET 1 

231  . 

Nl  x  Ml  - | 

2  32. 

CALL  REH(R(| ,M| .M.N.O.PTSETI 

733. 

IF (Ml  .GT.  11  GO  TO  20 

234. 

RETURN 

235- 

END 

236. 

ccc 

ccccccc 

237. 

c 

73B. 

c 

REM(R(L. mi .M.N.O. PTSET IS-  MAKES  a  HEAP  of  ELEMFNTS 

of  PTSFT 

239. 

r 

STA9TJNG  FROM  L  T 0  Ml.  PTSET  CONTAINS  M  POINTS 

IN  N 

240. 

c 

DIMFNSIONAL  SPACF. 
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2A|  . 

c 

2A2. 

CCCCCCCCCC 

2A3. 

SUBROUTINE  PEH1P(L. 

244  , 

Integer  »tseti 110.6 

245. 

integer  elrg  .x 

246. 

1  *  L 

247. 

J  *  2»  1 

2*8. 

X  a  Ml  1 

249. 

33  5  LL *  1  .N 

2S0. 

5 

R TSET (  x  ,LL  1  *  PTSET 

2S|  . 

EL4G  *  1 

252. 

lcIJ  .GT.MII  RETURN 

253. 

40 

I E { J  . GE  «  Mil  GO  TO 

254. 

C4LL  C0MP4RI J, J*l  , N 

255. 

IE  II  COMP  .UT.OI  J  1 

25f  . 

10 

C4LL  C  0M»  4PIX.J.N.0 

257. 

1=  II  COMP  .ge.o*  GO 

258. 

03  I  00  LL  *  1  ,N 

259. 

1  00 

PTSET  I  1  ,l_L  1  «  PTSET 

260. 

I  =  J 

26|  . 

J  *  24  1 

262. 

03  209  LL  *  1 ,N 

263. 

200 

PTSET!  |  ,.u  >  PTSET  1 

264  . 

GO  TO  30 

265. 

20 

flag  *  0 

266. 

30 

IE  (J  ,LE.  Ml  .AND. 

267. 

RETURN 

268. 

end 

269. 

CCCCCCCCCC 

270. 

C 

27  1  . 

C 

C0MP4RI I . J.N.O.PTSET . ] 

’2, 
3. 
<74. 
2  75. 

276. 

277. 

278. 
27q. 
280. 
281  . 
282. 
283. 
284  • 

285. 

286. 
287  . 
288. 
269. 
2»0. 
291  . 
2«2. 

293. 

294. 

295. 
’06. 
297. 

296. 

299. 

300. 


1  0 

0.  PT! 
>  J*  1 


r0  20 


FLAG  .50.11  GO  TO  40 


.  _  _ *Pi:-  COMP 4 RES  TWO  POINTS 

C  4  HANDLE  VECTOR  IN  0.  AND  RETURNS  IN  T COMP : 

C  1  IE  !TH  POINT  >  JTH  POINT, 

C  0  IE  ITM  POINT  «  JTH  POINT. 

C  -(  IE  ITh  POINT  *  JTH  POINT. 

CCCCCCCCCC 

SUBROUTINE  COMP  API  I . j.n.O.PTSET. ICOMPI 
INTEGER  “TSETI  I  IO.61.OI6t 
II  *  I 

100  III  *  01  l  I  I 

IE  (OTSETI I ,1 11  *  ,NE.  PTSETI J.I I  I  * »  GO  TO  10 
IE  (II  .LT.  Nl  GO  TO  5 
I  COMP  =  o 

peturn 

5  II  «  I l ♦ I 

GO  TO  100 

10  !E(PTSETI  I  .  1  I  I  »  .lT.  PT  SET IJ.1TI  II  ICOMP  z  -\ 

IE  IPTSETI  I  .1  I  I  I.GT.PTSETIJ. I  II  I »  ICOMP  «  I 

RETURN 

END 

CCCCCCCCCC 

c 

r  SXCHANt  I  .J.N.PTSET IS-SWITCHCS  the  positions  OF 
C  J  IN  ARRAY  PT  SET  WHICH  contains  A  POINT  SET 

C  SPACE, 

c 

CCCCCCCCCC 

SUBROUTINE  EKCHANI I . J.N.PTSETl 

integer  ptsfti i io.oi 


i  t 


W.R.T 


TwO  POINTS  I.  D 
IN  N  01  MENS 1 ON 4L 
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301  . 
10?. 
303. 
304  . 

305. 

306. 
30?  . 
30B  • 
30®. 

319. 

311. 

312. 

313. 

314. 

315. 

316. 
31?. 
318. 
310. 

320  . 

321  . 
3??. 
32  3. 
3?4. 

325. 

326. 
3??. 

329. 

320. 
?30. 

31  . 
j32. 
333. 
334  . 

335. 

336. 

337. 
33(1. 

330. 
340. 
34  1. 

342. 

343. 

344. 

345. 

346. 

347. 

348. 
340. 
350. 
351  . 

352. 

353. 
354  , 

355. 

356. 

357. 
359. 

359. 

360. 


03  100  II  =  1  .N 
I  TEMP  =  PTSET ( I  .1  I  I 
OTSETM.in  x  PTSCMJ.III 
RTSCTIJ.  Ill  *  ITEM  P 
100  CONTINUE 
RETURN 
ENJ 

CCCCCCCCCC 

c 

C  ENCODE ( MM. M. N.O .PTSET , NUM , I FL I t -GENERATES  GMMM  COOES  FOP  4 
C  SET  PTSET.  CONTAINING  M  POINTS  IN  N  0I»*ENSI0N4L  SPACE. 

C  NUM  WILL  CONTAIN  THE  NUMBER  OF  SUCH  Ghmm  ENCOOEO  TU»LES  IN  IT. 

C  TEL  IS  A  FLAG  TO  TUPN  OFF  ANO  ON  THE  POINTING  OF  INTERMEDIATE 

C  RESULTS. 

C  CALLS! 

C  (II  MtNMAX:  TO  PINO  THE  MINIMUM  C  MAXIMUM  EL  F MFN?  IN  A 

C  ’ARTICULAR  COLUMN  FOR  SOME  POINTS  CONSIDERED  TO  BE  IN 

C  CURRENT  HYOERCUBE. 

c 

CCCCCCCCCC 

SUBROUTINE  ENCODE (Mn.m, N. O. PTSET . NUM. I FL  » 

INTEGER  ’ ?  SE  T (  1  10.6|.0(6» .CODE! 1 2  I 
I F (  IFL  .EO.  II  GOTO  44 
PRINT  I 

I  FORM  AT ( // I 

PRINT  5. MM 

5  rORMAT(|ix.*GH *.!?.•  CODES  ARE:* I 

44  IF  (MM  .EO.  I  I  GO  To  10 

1=(MM  .EO.  N*ll  GO  TO  20 
N'JM  *  0 
N1  x  MM— | 

«»1  =  M  ♦  I 

I  =  I 

J  =  1 

DO  100  II  =  2. Ml 

IF  (II  .GT.  Ml  GO  TO  40 
FLAG  =  0 

30  200  31  *  1  .N| 

JJ  =  0(JII 

I r (PTSCT (  1 1 -( . JJ  |  ,EO.  PTSET I | I  .JJ! I  GO  To  ?00 
F.AG  x  t 
GO  TO  30 

200  continue 

30  IF  (FLAG  .EO.  1  I  GO  TO  *0 

J  x  j  ♦  | 

GO  T  3  100 

40  DO  300  Jl  x  I . N1 

J  J  x  0(  Jl  I 

C  DDE ( J 1  I  *,  PT  SET (  I , J J I 
300  CUNTINUE 

L  x  Nt 
NN  x  N  -  N1 
03  430  Jl  x  MM.N 
J  J  =  01  Jl  I 

CALL  MINMAXl’TSET.M.JJ,  I  .J. MIN. MAX  I 
L  =  L  ♦  I 
CODE  1 L I  x  MIN 
AOO  C  ODE  1 L  ANN  I  x  MAX 

N2  x  L 
N3  *  L  ♦  NN 
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161  . 
;62. 
36’. 
36*  . 

365. 

366. 
367  . 

366. 
360. 
370. 
37|  . 

372. 

373. 
37*. 
3  75. 

376. 

377. 

378. 
3  70. 
330. 
351  . 

382. 

383. 
38*. 

385. 

386. 

367. 
388  . 

389. 

390. 
”  I  . 

2. 

-<3. 

39*. 

395. 

396. 

397. 

398. 

399. 
*00. 
*01  . 
*02. 
*03. 
•  0*  . 
*05  . 
*06. 
*07  . 
*08  . 
*09. 
*10. 

•  It. 
*12. 
*13. 
*1*. 
*1  5. 
*15. 

*  I  7. 
*18. 
*19. 
*20. 


! 
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|F( |FU  ,F0.  t I  GO  TO  85 


PRINT  |  1  . 1C00EII  1  .  I  =  1  .  N1  » 

1 1 

FOP* AT <*0*8 

1  OX . • ( • . 3 1 3  1 

PRINT  | ’ ,< C30E( I # . I  =MM.N2I 

13 

FflPMAT  <  •  9 •  ♦ 

N ??  *  N?  A  1 

2  2  X , • ;*,313» 

PRINT  15. (CODE (  I  1 .1*  N22.N3I 

1  5 

FDPM  AT  i  •  A  •  • 
PKINT  |6 

33X. • ; 31 3» 

l  6 

FQPM  AT  (  '  A  •  • 

44  X •  •  |  •  I 

85 

I  *  !  1 

J  =  II 

MUM  *  NUN  A 

100 

C0^T 1MUE 
RE7URV 

1  0 

DO  1  I  1  I  *  1 

•  N 

JJ  =  0(11 

C  *-  L  M)NKUI  PT  SETT  .M.  JJ,  I  , M.M| N.MAXI 

cooe (ii  =  m (n 

coo? (  1 «mi  *  m * x 

l  t  I  CriNT  1  NUE 

|F(|FL  .EO.  I  I  RETURN 
JOINT  3. ( CODE!  !  I  .1 =1  INI 
8  FORM *T ( *0 • . 10 X . • ( • , 4 1 3» 

NN  =  N  »  1 

N2  =  2*N 

point  90 . ( CODE!  I  * • I  =  NN. M2 » 
oo  =3#«<TI  '**,  2  4X  # •  ;a.4!3l 

POINT  93 

03  F[IR«IT(  »♦  •  .  38X,*  >•> 

netuon 

20  0”  2?2  I  *  I  .  *0 

on  220  J  *  I  .  N 
JJ  «  OIJI 

220  CODEIJI  =  PTSETM.JJI 

|C  (  I  El  .EO.  I  I  GO  TO  222 
JOINT  22. (CODE! J 1 . J  =  I ,N» 

22  F  O  pm *Tf  *0*  ,IOX,,(#,*t3l 
point  23 

23  E  nPM *  T (  • ♦ • ,  2*X , *  »•» 

222  CONTINUE 

RETURN 
END 

CCCCCCCCCC 

c 

r  MI NMA* ( PT SET , M . | COL . I . J. MIN.MAX I : -  TO  FIND  THE  MINIMUM  t 

C  MAXIMUM  ELEMENTS  OE  COLUMN  ICOL . ST*PTING  EP»M  I  TO  J  IN  *PR *V 

C  PISET. 

c 

CCCCCCCCCC 

SU8POUT INF  M|NMAX(PTSET,M, ICOL. I.J. MIN. MAXI 

integer  ptseti i 10.61 

MIN  =  PT  SE  T  <  I . I  COL  I 
MAX  S  PTSFT I  I .  lent  I 
33  100  II  =  I . J 

I  1  (’TSETIII.ICOLI  ill*  MINI  M|N  =  RTSET I  II  .  ICOL I 
tr  IOTSET (  I  I  .  ICOL *  .GT.  maxi  max  *  PTSET  I  I  I  .  ICOL  I 
too  CONTINUE 

return 

ENO 
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h 

i 


'1 

i 


I 


♦  21  . 
♦  22. 

♦  23. 

♦  ?♦. 

♦  25. 

♦  26. 

♦  2T. 

♦  28. 

♦  29. 

♦  30. 

♦  31  . 

♦  32. 

♦  33. 

♦  34. 

♦  35. 

♦  3ft. 

♦  37. 

♦  38. 

♦  39. 

♦  ♦0. 
♦  ♦1  . 
♦  ♦2. 

♦  ♦3. 

♦  ♦♦. 

♦  ♦5. 

♦  ♦6. 

♦  ♦7. 
♦♦8. 

♦  ♦9. 

♦  50  . 

51  . 

2  . 

»S3. 

♦  54. 

♦  55. 

♦  56. 

♦  57. 

♦  58. 

♦  59. 

♦  60. 
♦  61  . 
♦  62. 
♦63. 

♦  64. 

♦  65. 

♦  66. 

♦  67. 

♦  68. 

♦  69. 

♦  70. 

♦  71  . 

♦  72. 

♦  73. 

♦  74. 

♦  75. 

♦  76  . 

♦  77. 
478. 
4  79. 

♦  80. 


CCCCCCCCCC 

c 

C  HEURIS(sTSET.M.N.NCODE. IFL1 :-  TO  STORE  THE  NUMBER  OF  GH  t  FNCODED 

C  TUP.FS  FOR  1  =  1.2 . NAI  IN  ARR4Y  NCOOE.  THE  SET  OF  HANDLES 

C  IS  CHOOSEN  OY  HEYRISTIC  METHOD. 

c  calls: 

C  (II  PBOCES:  TO  COMPUTE  the  measure  of  priority  P|I»  FOR  EACH 

C  COLUMN  I. 

c  <2i  sort:  to  get  the  proper  set  of  handles  corresponding  to 

C  to  BFST  VALUES  of  their  PRIORITY  P(t». 

C  (31  SORTPT:  TO  SORT  POINT  set  ptset  K.R.T.  CHOOSEN  SET  OF 

c  handles. 

C  (41  encode:  TO  GENERATE  GH  ENCODED  TUPLES, 

c 

CCCCCCCCCC 

SUBROUTINE  HEURIS(PTSET.M. n.ncode. IFL 1 
INTEGER  PTSET (  110.6  1.0(6  I.PI61.NC0DE (5. 31 
CALL  PROCESIPTSET.P.M.N! 

PRINT  1 

1  F  OPMA  T ( / / / *  I  1 X,  1  PR  INTOUT  OF  ARRAY  P  IS:*1 
PRINT  2  . ( P11 I  .1 =1 .Nt 

2  rORMAT ( 10X .51 1 0 1 
DO  200  I  =  I • N 

200  0111  =  I 

CALL  SORT ( P .0 . N I 
PRINT  3 

3  =OPMAT ( «0*  . tox . ’HANDLES  CHOSEN  ARE  5 • » 

PRINT  4  ,(0( I  1 .  1  =  1  ,N> 

♦  FJRMATM  *  .  1  OX  .4151 
CALL  SORTPT (PTSET .O.M.Nt 
I F ( I=L  .EO.  II  GO  TO  ♦♦ 

PRINT  5 

S  format (* 0* . 10X ,• SORTED  POINT  SET  IS:*1 

DO  300  I  =  I .M 

300  PRINT  I  0 • ( PTSET (I . J I . J=  l.NI 

10  FORMAT (|  ix  .* ( •  .413.*  »•  » 

PR  I NT  II 

11  FORMAT(/7/.20X. ‘FOLLOWING  ARE  GH  CODES  V.R.T. ABOVE  HANDLES:* » 

♦  ♦  Nl  *  N  ♦  I 

DD  400  1  =  l.NI 

CALL  ENCODE! I . M . N, 0. PTSET, NUM. IFL! 

1=1  1  .LB • I  .OR.  I  ,GE.  Nil  GO  TO  400 

NCOOE (  I  —  I , I  1  =  NUM 
NCOOE ( 1-1.21  =  NUM 
NC30EI I- I .31  =  NUM 
♦00  CONTINUE 

RETURN 
END 

CCCCCCCCCC 

C  EXHAJS(PTSET. M.N, NCOOE. IFLi:-  TO  STORE  the  MINIMUM  NUMBER  OF  GH I 

C  ENCODED  TUPLES  FOR  l=l.2,...NM  IN  ARRAY  NCOOE. ALL  THE 

C  POSSIB.E  SET  OF  HANDLES  ARE  USED  TO  GENERATE  GH  CODES  t  THE 

C  ONE  WITH  MINIMUM  NUMBER  OF  TUPLE  IS  TAKEN  AS  OPTIMAL  GH 

C  CODES, 

c  calls: 

r  in  pcdmb:  to  fin  all  possible  sets  of  handles. 

C  (21  S3  2  TP  T :  TO  SORT  THE  P01  NTSFT  W.R.T.  CURRENT  SET  OF  HANDLES 

C  131  FNCODE:  TO  GENERATE  GH  ENCODED  TUPLES. 

C 


f 


I 

) 


J 

I 

\ 

i 


! 

i 

t 

( 
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y 


481. 

20 

RETURN 

482. 

10 

DO  222  1  =  l.N 

483. 

JJ  *  0(  I  ) 

484. 

CALL  MI NMAXIPTSET. ORDER. JJ.  1 .NPTS .M I N.MAX 1 

48S. 

CODE ( I )  =  MIN 

486. 

CODE! I4N)  =  MAX 

487. 

222 

CONT 1NUE 

488. 

DEN  =  DENSTY (CODE. MM .N.l  .NPTS) 

489. 

IE (DEN  .GE.  THRESH}  GOTO  444 

490. 

NUM  =  0 

491  . 

return 

492. 

444 

IF  1FLG  .EO.  1)  GO  TO  666 

493. 

PRINT  2 1 . (CODEC  I ) . 1  =  1 .N> 

494. 

21 

FORMAT 1 *0*. 10X. •( >.613) 

495. 

N 1  =  N+ 1 

496. 

N2  =  N4  N 

497. 

PRINT  22. (CODE! !> . I=N1,N2> 

498. 

22 

FORMAT! *4' ,30X. • I *.613) 

499. 

PRINT  23 

500. 

23 

FORMAT ( *4*.S0X. • } • ) 

501. 

666 

NUM  =  1 

502. 

RETURN 

503. 

END 

504. 

CCCCCCCCC 

505. 

c 

506. 

C  I 

MINMAX(PTSET  .ORDER.  J  J.  I.  J.  MIN. MAX}  FIND  THE  MIMMUK  C  MAXIMUM 

507. 

c 

OF  ELEMENTS  OF  COLUMN  JJ  .STORED  IN  PTSET. THEIR  POINTERS  ARE 

508. 

c 

STORED  In  ARRAY  OROER  STARTING  FROM  I  TO  J. 

509. 

c 

\|0. 

cccccccccc 

11  . 

SUBROUT INE  MINMAXfPTSET. ORDER. JJ. I. J. MIN. MAX) 

12. 

INTEGER  PTSET (1 1 0 .6} . ORDER ( 110 ) 

13. 

1 I  =  ORDER (  I  I 

514. 

MIN  *  PTSET! I I. JJ} 

SIS. 

MAX  =  PTSET ( T I .JJ} 

516. 

JO  130  K  =  t.J 

517. 

1 I  =  ORDER (K1 

518. 

IF(PTSETI I I . JJ)  .LT.  MIN)  MIN  =  PTSET(II.JJ) 

519. 

IF!  PTSET (11 • J J  )  .GT.  MAX)  MAX*  PTSET(Il.JJ) 

520. 

100 

CONTI NUE 

521  . 

RETURN 

522. 

END 

523. 

cccccccccc 

524. 

c 

525. 

c 

OENSTY(COOE<MM, I . J) J-  COMPUTES  THE  DENSITY  OF  HYPERCUBE 

526. 

c 

GENERATED  BV  GHM  ENCODED  TUPLE  STORED  IN  CODE  WHICH  INCLUDES 

527. 

c 

points  pointed  to  by  pointers  in  order  starting  from  i  to  j. 

528. 

c 

529. 

cccccccccc 

530. 

FUNCTION  DENSTY(CODE.MM.N.I.J) 

531  . 

INTEGER  CODE! 121 

532. 

NUMPTS  *  J- 1  41 

533. 

VOL  *  1 

534. 

Ml  *  MM  4  1 

535. 

JO  100  X  *  Ml  .  N 

539. 

VOL  =  VOJ.MCODE!  K4N-MM) -CODE  !K  >41  ) 

537. 

100 

CONTINUE 

538. 

DENSTY  *  NUMPTS  /  VOL 

539. 

RETURN 

540  . 

END 
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54  1  . 

542. 

543. 

544. 
54  5. 

546. 

547. 
54  A. 

549. 

550. 
551  . 

552. 

553. 

554. 

555. 

556. 

557. 

558. 

559. 

560. 
561  . 
562. 

3. 

4. 
565. 
566  . 

567. 

568. 

569. 

570. 
571  . 

572. 

573. 


L  *  L  ♦  t  TERM 
I  I  *  I!  ♦  1 
LI  »  L  -  1 

IP  1  I  1  .GT.  LI  I  GO  TO  20 
00  500  1 2  -  I  I  .  Lt 
T  ABLE (12*  COL  I  «  T  ABLE (If. COL  I 
500  CONTINUE 

20  IP  (I  .GT.Nlt  SO  TO  400 

CALL  FILLdABLE.COL.N.  1I.LI 
400  CONTINUE 

300  CONTINUE 

K  *  NUN  *  K 
NUN  *  NUN  -  1 
<"OL  =  COL  ♦  1 
IP  (NUN  .GE.2I  GO  TO  10 
RETURN 
ENO 

CCCCCCCCCC 

c 

C  P I L. ( TABLE . COL. N. If. LI**  TO  STORE  SONE  ELEMENTS  IN  LTM  ROW  OP 

C  ARRAY  TABLE. 

C 

CCCCCCCCCC 

SUBROUTINE  P 1 LL I T ABLE . COL. N. II .LI 
INTEGER  t ABLE (720.61 
INTEGER  COL.  COL  I 
COL  I  *  COL  ♦  I 
I  X  =  TABLE! 1 1 .COL  I 
00  100  1  *  COL  I .N 

100  T  ABLE  (L  .  1*1  I  *  TABLEdI.ll 

TABLEfL.NI  *  IX 
RETURN 
ENO 


APPENDIX  B 


1.  ✓•JOB PAR*  T=I2.0).L=4 

2.  //  EXEC  FORTGCLG 

3.  //FORT.SYSIN  00  • 

4.  CCCCCCCCCCC 


5.  C 

6.  C  MAIN  PROGRAM  TO  OPTIMIZE  GH  ENCODING  FOR  A  PCINTSET. 

7.  C  SUCH  THAT  EACH  HTPERCUSE  HAS  DENSITY  NO  LESS  THAN  A  GIVEN 

8.  c  threshold  density. 

9.  c  ptset;  array  containing  pointset. 

to.  c  thresh:  threshold  density. 

it.  c  m:  no. of  points  in  point  set. 

12.  c  n:  dimensional  space. 

13.  c  flg:  flag  used  TO  TURN  off  c  on  printing  of  intermediate 

14.  C  RESULTS. 

is.  c  calls 

16.  c  (1)  GENPTS:  TO  GENERATE  M  POINTS  IN  N  DIMENSIONAL  SPACE  t 

17.  C  STORE  THEM  IN  ARRAY  PTSET. 

18.  C  (2)  optmum:  to  GENERATE  OPTIMAL  GH  ENCODED  TUPLES.  AND 

19.  C  RETURNS  THEIR  NUM3ER  IN  NCODE.  IF  FLG  IS  ON  .PRINTS 

20.  c  those  gh  encoded  tuples  also,  thersh  is  the  threshold 

21.  C  DENSITY  SUCH  THAT  EACH  HYPERCUBE  HAS  DENSITY  NO  LESS 

22.  C  THAN  THRESH. 

23.  C 

24.  CCCCCCCCCCC 

25.  INTEGER  PTSET! 1 10.61 .TABLE ( 1 1) 

26.  INTEGER  FLG 

2T.  N  =  3 

28.  Fi_G  =  I 


1 . 

DJ  300  M  =  10.103.10 

.2. 

CALL  GENPTSCPTSET.M.NI 

32.  1 

Du  400  ITH  =  1 0. 1 10. 10 

32.2 

NJEN  =  ITH  -  10 

32.3 

THRFSH  *  NDEN/100.0 

33. 

IF  (FLG  .EO.  1 )  GO  TO  10 

34. 

PRINT  1 

35. 

1 

F JRMAT ( •  1  • ,  10X  .  1  INITlAIL  POINT  SET  is: 

36. 

Ou  100  I  =  1 .M 

37. 

PRINT  2  .(PTSET! I .J> ,J»1 ,N> 

38. 

o 

o 

CUNTINUE 

39. 

2 

FuRMAT ('O'  .I0X.1!1  .613.*  )'  ) 

40. 

PRINT  3 

41  • 

3 

FURMAT (////. 100! I  > 

42. 

PRINT  4 

43. 

4 

FURMAT (2/.25X, • HEUR I  ST IC  APPROACH'  I 

44. 

PRINT  5 

45. 

5 

F  JPMAT  t  25X  *  • - - »  ) 

46. 

1  0 

CALL  OPTMUM(PTSET.M.N.THRESH.NCOOE.FLG) 

47. 

400 

TAOLEI  I  TH/101  =  NCODE 

49. 

PRINT  6,N 

50. 

6 

FURMAT(///,10X, 'DIMENSION  =  *.I5) 

51  . 

print  t.m 

52. 

7 

FuRMAT  ( 10X  ,' NO  .  OF  POINTS  IN  SET=  '.151 

53. 

PRINT  b 

54  . 

0 

FORMAT (• 0 ', 10X ,' threshold  ',10X,'»  of  c 

55. 

DJ  20  0  I  =  I  • 1  l 

56. 

ITH  =1*10-10 

57. 

PPINT  9.  1TH.TADLEI  I  1 

58. 

200 

CUNTINUE 

59. 

300 

CJNTINUE 

60. 

9 

FuRMAT (t0X.I10.I0X.ll0) 

ODES  GENERATED' > 


23 
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STOP 

END 

ccccccccccc 

c 

C  RAN  DOM (  X )  :  GENERATE  A  RANDOM  NUMBER  X .  J  <=  X  <=1.U.  WITH 

C  UNIFORM  U IS TRi dUT I  ON. 

C 

CCCCCCCCCCC 

SUBROUTINE  RANDOM! X) 

R  £  AL  X 

INTEGER  A/1D727/.  MULT/2S211/.  E ASE/U2766/ 

A  -  MODI MULTWA, BASE J 

X  =  ECOrT I  AJ/FCQAT I  BASE  I 

RETURN 

END 

cccccccccc 

c 

C  GtNPTSI PT S ET , M, N) : -  TO  GENERATE  TEST  DATA  SET  IN  ARRAY 

C  P1SET  i  WHICH  CONTAINS  M  NUMBERS  IN  N  DIMENSIONAL 

C  SPACE. 

c  calls: 

c  (i)  randum:  tc  get  a  random  number. 

c 

cccccccccc 

SUBROUTINE  GtNPT SI PTSET ,  M.N) 

INTEGER  PT  SET (110.61 
DO  100  J  =  1.  N 
DO  ICO  1  =  1  .  M 

CALL  RANOUMIU) 

PTbET(I.J)  =  I  NT  1 9  *  U»  +  1 

ioo  continue 

RETURN 


END 

cccccccccc 

c 

OPT MUM ( Pis  ET  , M. N,  T  HR"Sm, MCuDE ,FLG  > 

TO  GENERATE  OPTIMAL  SH  LNCuO  1  NG  FOR  A  oIVEN  THRESHULO 

density,  thresh;  and  retuans  the  number  of  such  tuples 

I  N  NC  ODE  • 


calls ; 
1 1  > 


c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

cccccccccc 

SUuKOUT  INC  DPT  MJ  Ml PTSET .M.N  .  T  UR £ SH . NEC DE  .FEE ) 

INTEGER  P1aET(110.6J.P|6J.Q(61 .ORDER! 1 10 ) 

intecer  flg 

NCODE  =  0 

DC  100  1  =  l .M 


(  2) 
(  o> 
(  A  ) 


PHCCEb:  TO  COMPUTE  THE  MEASURE  DF  PRIORITY  PI  I)  FOR 
EACH  COLUMN  I . 

SGPT:  TO  GET  THE 

TO  BEST  VALUES 
SORTPT:  TO  SORT 

HANOLES. 

encode:  tg  generate  ghm  encoded  tuples  with  density 

no  LEbS  than  threshold  DENSITY. 


PROPER  SET  OF  HANDLES  CORRESPONDING 

of  their  priorities,  phi. 

POINTSET  IN  PTSET  a.R.T.  NE»  SET  UF 


100 

(11=1 

NPTS  -  M 

KM  =  l 

1  1  0 

JFL  =  0 

DO 


}00  I  =  1  .N 
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1  21  . 
122. 

123. 

124. 

125. 

126. 

127. 

128. 

129. 

130. 
131  . 

132. 

133. 
1  34. 
1  35. 

136. 

137. 

138. 

139. 

140. 
141  . 

142. 

143. 

144. 
1  45. 

146. 

147. 
1  48. 

9. 
,0. 
1S1  . 
152. 
I  S3. 
1S4. 


157. 

158. 

159. 

160. 
161. 
162. 

163. 

164. 

165. 

166. 

167. 

168. 

169. 

170. 
171  . 

172. 

173. 

174. 
1  75. 
1  76. 

177. 

178. 

179. 

180. 


500  0(1)  =  1 

IF  ( MM  .GT .  1 )  GOTO  80 

GALL  ENCODE (MM. N. PT SET. ORDER . 0 . NPTS . 1 FL .NUM .THRESH , F LG ) 
IF  (NUM  .EO.  0)  GOTO  70 
NCOOE  =  NUM 
RETURN 

7  0  M  M  =  MM  +  I 

80  CALL  PROCE  S(PTS£T . ORDER. P«  NPTS . N) 

IF  (FLG  .EO.  I )  GO  TO  30 
PRINT  l.(P< 1 ) .1=1 .N) 

1  FORMAT!//. 10X. 'PRINTOUT  OF  P  ARRAY  ISl*.  615) 

30  CALL  SORT!  P.O.N) 

IF  (FLG  .EO .  1 1  GO  TO  40 

PRINT  2  .  (Q(  I).  1  =  1 .N) 

2  FORMAT ( // • 1  OX . ' ORDER  OF  VARIOUS  HANLES  IS  5  •  .  615) 

AO  CALL  SURTPT (PTSET.O. ORDER. NPTS, N.M) 

CALL  ENCODE (MM. N. PT SET .ORDER. 0. NPTS. 1 FL .NUM .THRESH. FLG ) 
NP  =  0 

DO  200  1  =  l.NPTS 

1 F ( ORDER ( I )  .EO.  0)  GO  TO  200 

NP  =  NP  +  1 

ORDER (NP)  =  ORDER ( I ) 

200  CONTINUE 

NPTS  =  NP 

NCOOE  =  NCODE  ♦  NUM 
IF  (NPTS  .EO.  1)  GO  TO  20 
IF  (NPTS  .LE.O)  RETURN 

IF  (NUM  .LT.l  .AND.  I FL  .EO.  I  )  MM  =  MM  ♦  1 

IF  (NUM  .LT.  1  .AND.  I FL  .EO.  0)  GO  TO  20 
GO  TO  110 

20  IF  (FLG  .EO.  1)  GO  TO  SO 

PRINT  5.  MM 

5  FORMAT  (// ,1  OX HANDLE  LENGTH  =  ' . 1 5 ) 

PRINT  6  . ( 0{  11.1=  I  . MM) 

6  format (✓ .iox . 'handles  are  :*.6is) 

PRINT  7 

7  FORMAT!/. 1  OX . •  CORRESPONDING  GH  CODES  ARE 1 • ) 

50  NCOOE  =  NCOOE  ♦  NPTS 

IF  (FLG  .EO.  II  RETURN 
DO  300  I  =  l.NPTS 
I  I  =  ORDER! I  1 
00  333  J=I.N 
J J  =  0( J) 

333  P(J)  =  PTSETUI.JJI 

POINT  8. (P( K) .K= I , MM  ) 

8  FORMAT ( 1  OX .'  (  '  .615. •  )•  ) 

IF  (N-MM  .EO .01  GO  TO  300 

N I  =  MM+1 

PRINT  9  .(P(  K)  .K=M  ,N) 

9  FORMAT!  ,42X,  • J  • ,51  5) 

PRINT  |0  ,(P(  K 1 .K=NI  ,  M 

10  FORMAT! ,65X.' I*. 51 51 
PRINT  II 

11  FORM  AT (  • ♦»  ,85X, •  I •  I 

300  continue 

RETURN 

ENO 

cccccccccc 

c 

C  PROCc.  S( P  TSET , OROE  R.NPTS.N! I—  COMPUTES  THE  MEASURE  OF  PR1CRTY 
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C  ,3(1)  FOR  EACH  column  1  WHERE  1=1... ..N.  FOR  A  GIVEN 

C  POINT  SET  IN  ARRAY  PTSET  CONTAINING  NPTS  NO.  OF  POINTS 

C  IN  N  DIMENSIONAL  SPACE. 

c  calls: 

c  (i)  hpoder:  to  compute  p( i )  for  one  particular  column  i. 
c 

cccccccccc 

SUdROUT 1NE  PROCES( PTSET . ORDER ,P , NPT S . N > 

IN  TEGER  PTSET! 1  10.6) .ORDER!  110) . TEMPR! 1 10I.PI6I 
00  100  I  =  l.N 
DO  200  J  =  l.NPTS 
JJ  =  ORDER (J ) 

200  TEMPR(J)  =  PTSETIJJ.il 

CAcL  HPODER!  I  .TEMPR, NPTS.P) 

IOC  CONTINUE 
RETURN 
END 

CCCCCCCCCC 

c 

C  HPODER!  I.  TEMPR, M. PH-  COMPUTES  P(I)  FOR  A  COLUMN  VECTOR  TEMPR 
C  .IT  USES  HEAP  STRUCTURE  FOR  THAT  AND  HENCE. 

C  CALLS: 

c  !i)  reheap:  to  make  heap  of  elements  of  column  vector  tempr. 
c 

cccccccccc 

subroutine  hpoder ( i .tempr. m.p > 

INTEGER  TEMPR ( 1 1 0 ) ,P (6 ) 

integer  temp. elment.psum, sum 

Ml  =  M 
L  =  M/2  ♦  l 
10  L  =  L-l 

CALL  REHEAP(L .Ml .TEMPR ) 

IF  !-  .GT.  I)  GOTO  10 
ELMENT  =  0 
P  SUM  x  0 
SUM  =  0 

20  TEMP  =  TEMPRI1 ) 

TEMPR! 1)  =  TEMPRIM1 1 

IF  (CLMENT  .EO.  TEMP»GO  TO  30 
IFISUM  .GT.  1  »  PSUM  =  PSUM  ♦  SUM  *  SUM 
CLMENT  =  TEMP 
iU«  =  1 
„0  TO  40 

30  SUM  =  SUM  -»■  1 

40  Ml  =  Ml  -  1 

1FIM1  .GT.  0>  GO  TO  SO 

IF  (SUM  .GT.  II  PSUM  =  PSUM  ♦  SUM  »SUM 

P( I )  =  PSUM 

RETURN 

50  CALL  REHEAP! I. Ml  TEMPR) 

GUTO  20 
END 

cccccccccc 

C  REHE AP!L . M l , TCM=R ) :-  MAKES  A  HEAP  OF  THE  ELEMENTS  OF  ARRAY  TEMPR 
C  STARTING  FROM  T E “PR ! L )  TO  TEMPR(Ml). 

C 

cccccccccc 

SJoROUTINC  PEHE AP ( L. Ml  ,  TEMPR ) 

1 N  f  E G£  R  TEMPR ( l  I  0  I 

integep  flag  .  x 


ft/ 


/ 


2*1  . 
2*2. 
2*3. 
2*«. 
2*5. 
2*6. 
2*7. 
2*6. 
2*9. 

250. 

251. 

252. 

253. 

254. 

255. 

256. 

257. 
25B. 

259. 

260. 
261  . 
262. 
263. 
26*. 

265. 

266. 

267. 

268. 
'9. 

0. 
271  . 

272. 

273. 
27*. 

275. 

276. 

277. 
276. 

279. 

280. 
281  . 
282. 

283. 

284. 

285. 

286. 

287. 

288. 

289. 

290. 
291  . 

292. 

293. 
29*. 

295. 

296. 

297. 

298. 

299. 

300. 
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1 1  =  L 
JJ  *  2*11 
X  =  TEMPR (111 
PL  AG  -  1 

IF  (JJ  .GT.  Ml)  RETURN 
*0  1F1JJ  .GE.  Ml)  GOTO  10 

IF  (TEMPR(JJ)  .LT.  TEMPR ( J J  +  1  )  )  JJ  =JJ+1 
10  I F 1 X  .GE.  TEMPR(JJ))  GOTO  20 

TEMPR (II)  *  TFMPR(JJ) 

I  (  *  J  J 
JJ  *  2*11 
TEMPR ( 1 1 )  *  X 
GOTO  30 

20  FLAG  =  0 

30  IF  (JJ  .LE.  Ml  .AND.  FLAG  .EO.  1)  GOTO  *0 

RETURN 
ENU 

cccccccccc 

c 

C  SORTIP.O.N >:-  SORTS  p  ARRAY  CONTAINING  N  ELEMENTS  AND  DURING 

c  this  sort  also  moves  the  elements  of  array  o  in  the  same  fash- 

C  ION  SO  THAT  INORDER  TO  CHOOSE  A  HANDLE  SET  OF  M  LENGTH.  ME  CAN 

C  CHOOSE  first  m-i  handles  stored  in  first  M-l  position  of  o. 

c 

cccccccccc 

SUdRUUT INE  SORTIP.O.N) 

INTEGER  0(6). p(6) 

INTEGER  FLAG.  TOP.  TEMP .  BOTTOM 
FLAG  =  1 

BOTTOM  «  N 
TOP  *  1 

10  1 F ( FLAG  .NE.  0  .AND.  TOP  .LT.  BOTTOM)  GU  TO  20 

RETURN 

20  FLAG  =  0 

I  I  =  BOTTOM  -1 
03  100  I  =  l.tl 

IF  !P(I)  .GE.  P(I*1)>  GO  TO  100 
X  *  P< 1  ) 

P( 1 )  =P1 1*1 1 
P( 1*1 )  =  X 
TEMP  =0(1) 

0(  I )  =0(  1*1  ) 

0(  1*1 )*TEMP 

flag  =  i 

100  CONTINUE 

BUTTON  =  BOTTOM  -1 

Gu  TO  10 

END 

CCCCCCCCCC 

c 

C  SDRTPTIPTSET. D. ORDER. M.N, M2 )1-  TO  SORT  POINT  IN  ARRAY  PT SET 

C  W.R.T.  NE*'  HANDLE  VECTOR  IN  O  .  ORDER  CUNtAINTS  POINTERS  TO 

C  POINTS  IN  PTSET  WHICH  HAVE  NOT  BEEN  DELETED  YET.M  IS  THE 

C  NUMBER  OF  SUCH  POINTS  WHILE  THERE  WERE  M2  POINTS  IN  THE 

C  ORIGINAL  POINT  SET.  N  IS  THE  DIMENSION  OF  SPACE. 

C 

C  HEAPSORT  IS  USED  HERE.  THEREFORE: 

c  calls: 

c  (  i)  rehip:  to  make  heap. 

c  (2)  compar:  to  compare  order  of  two  points. 
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301. 

c 

(3)  exchang:  to  exchang  position  of  two  points 

IN 

PTSET . 

302. 

c 

303. 

CCCCCCCCLC 

304. 

SU9K3UT  I  NE  SORTPT  (PTSE  T.O.  ORDE  R  .  M.  N  .  M2  ) 

305. 

INTEGER  PT  SET  1 1 1 0 . 6 1 . 0 1 6 1 . ORDER (  1 1 0 ) 

306. 

Ml  =  M 

307. 

L  =  Ml/2  ♦  1 

308. 

1  0 

L  =  L  -  1 

309. 

CALL  REHIP<L.Ml,M2 .N , 0 . PTSE T .ORDER ) 

31  0. 

IF IL  .GT .1 1  GOTO  10 

311. 

20 

CALL  E  XCH AN (ORDER! 1) .ORDER (Ml >.N. PTSET ) 

312. 

Ml  =  Ml-1 

31  3. 

CALL  REHIPd  .Ml  .M2  .N,  O.PTSET, ORDER  ) 

314. 

1FI Ml  .GT.  1 )  GO  TO  20 

315. 

RETURN 

31  6. 

END 

317. 

ccccccccc 

318. 

c 

319. 

c 

REHI P( L ,M I .M2 .N .0 , PTSET .ORDER) : -  MAKES  A  HEAP 

STRUCTURE 

320. 

c 

of  those  points  in  ptset  whose  pointer  are 

IN 

ARRAY  ORDER 

321  . 

c 

starting  from  order(l)  to  ordercmii. 

322. 

c 

M2  IS  the  NUMeER  OF  POINTS  IN  ORIGINAL  PCINTSET 

323. 

c 

324. 

cccccccccc 

325. 

SUdROuTINE  REH I P ( L • Ml  . M2 . N, 0. PTSE T . ORDER  1 

326. 

INTEGER  PTSET  < 1  I  0 . 6 1 . 0( 6 1 . ORDER ( 1 10  I 

327. 

INTEGER  FLAG  .X 

328. 

I  I  =  L 

29. 

JJ  *  24  I  1 

10. 

I  =  ORDER ( I 1  1 

11  . 

X  =  M2  ♦  1 

332. 

DO  5  LL=l  .  N 

333. 

5 

PTSET ( X .LL 1  -  PTSETCI.LL) 

334. 

FLAG  =  1 

335. 

1F1JJ  .GT.M11  RETURN 

336. 

40 

IF! JJ  .GE .  Ml  1  GO  TO  10 

337. 

CALL  COMPARinRDERt J J (. ORDER (JJ. 1 >.N. O.PTSET. 

ICOMP) 

338. 

IF  (ICO  MP  .LT.O)  JJ  =  JJ+l 

339. 

10 

call  C  OMP AR ( X.ORDER(JJ) .N .O . PTSET . I  COMP > 

340. 

IF  ( 1  COMP  • GE .01  GO  TO  20 

341  . 

I  =  ORDER ( I I 1 

342. 

J  =  ORDER ( J J 1 

343. 

DO  100  LL  =  l.N 

344. 

too 

PT  SET (  1 .LL 1  =  PTSET(J.LL) 

345. 

1  I  =  J  J 

346. 

JJ  =  24U 

347. 

1  =  ORDER (11) 

348. 

DO  200  LL  =  l.N 

349. 

200 

PTSET(I.LL)  =  PTSET ( X .LL I 

350. 

GO  TO  30 

351  . 

20 

FLAG  =  0 

352. 

30 

IF  (JJ  .LE.  Ml  .AND.  FLAG  . EO . 1  I  GO  TO  40 

353. 

RETURN  ' 

354. 

END 

355. 

CCCCCCCCCC 

356. 

c 

357. 

c 

COMPAR ( I , J .N. O.PTSET . ICOMP) J-  COMPARES  TWO  POINTS 

i 

1  J  W.R.T 

358. 

c 

handle  vector  in  o.  and  returns  in  icomp: 

359. 

c 

1  IF  ITH  POINT  >  JTH  POINT. 

360  . 

c 

v  IF  ITH  POINT  =  JTH  POINT. 
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361.  C  -1  IF  ITH  POINT  <  JTH  POINT. 

362.  C 

363.  CCCCCCCCcC 

364.  SUBROUTINE  COMPARI t . J . N . O.PT SET . I  COMP » 

365.  INTEGER  PTSET (1 10.61.0(6) 

366.  11  *  1 

367.  100  III  *  Dlltl 

368.  IF  IPTSET ( I . I  1 1  )  ,NE.  PTSET (J, I] I  ) )  GO  TO  10 

369.  IF  (II  .LT.  N)  GO  TO  5 

370.  ICUMP  =  0 

371.  RETURN 

372.  S  II  *  1141 

373.  GO  TO  100 

374.  10  1FIPTSETC 1.111 1  .LT.  PT SE T( J . 1 1  I  1 )  I  COMP  *  -1 

375.  IF  (PTSET (1.1II  i.GT .PTSET (J.1II  I)  I  COMP  =  1 

376.  RETURN 

377.  ENil 

376.  CCCCCCCCC 

379.  C 

380.  C  EX CHAN ( I . J.N. PTSET l SWITCHES  THE  POSITIONS  OF  TWO  POINTS  I  t  J 

381.  C  IN  ARRAY  PTSET  WHICH  CONTAINS  A  POINT  SET  IN  N  DIMENSIONAL 

382.  C  SPACE. 

363.  C 

364.  CCCCCCCCCC 

385.  SUBROUTINE  EXCHANl I . J.N. PTSET! 

386.  INTEGER  PTSET(110.6> 

387.  03  100  II  =  I . N 

366.  I TEMP  »  PTSET (1,1 I) 

389.  PT  SET ( I . I  I  I  *  PTSET ( J . I  I  I 

0.  PTSET(J.II)  <*  ITEMP 

*  I.  100  CONTINUE 

>2.  RETURN 

93.  END 

394.  CCCCCCCCCC 

395.  C 

396.  C  ENC0UE(M. N, PTSET. ORDER. O. NPTS. 1FL .NUM. THRESH. FLG) GENERATE 

397.  C  5HM  CGDES  FOR  A  POINT  SET  PTSET.  CONTAINING  NPTS  POINTS  IN 

396.  C  N  UIMENS1DNAL  SPACE.  ONLY  THOSE  HYPEMCUBES  WHICH  HAVE  OESITV 

399.  C  =>  THAN  THFESHOLO  OENS1TV,  THRESH.  ARE  INCLUDED  IN  OPTIMAL  GH 

400.  c  Encoding,  fop  a  ginen  m.  nun  contains  the  number  of  such 

401.  C  HYPERCUBE  generated  w.r.t.  a  GIVEN  handle  vector  IN  O. 

402.  C  ORDER  ARRAY  CONTAINS  POINTERS  TO  POINTS  IN  PTSET  WHICH  ARE 

403.  C  STILL  IN  POINTSET.  FLG  IS  A  FLAG  USED  TO  TURN  OFF  t  ON 

404.  C  PRINTING  OF  OPTIMAL  GH  ENCODED  TUPLES. 

405.  C  calls: 

406.  C  *  (11  MINMAX:  TO  FIND  THE  MINIMUM  L  MAXIMUM  ELEMENT  IN  A 

407.  C  P APT  1 CULAR  COLUMN  FOR  SOME  PONTS  CONSIDERED  TO  BE  IN 

408.  C  IN  THE  CURRENT  HYPERCUBE. 

409.  C  (21  DENSITY:  TO  COMPUTE  THE  DENSITY  OF  A  GENERATED  HVPCRCUBE. 

41 0.  C 

411.  CCCCCCCCCC 

412.  SUBROUTINE  ENC0DE(M, N, PTSET. ORDER. D. NPTS. IFL .NUM, THRESH. FLG I 

413.  INTEGER  FLG. FLAG 

*1*.  INTEGFR  P  TSET I  1  10.61.0(6) .C  ODE (121. ORDER ( 1 1  0  I 

415.  MM  x  m  — 1 

416.  IF  (FLG  .EQ.  11  GO  TO  70 

417.  PRINT  1 

418.  I  FORMAT (/✓  I 

419.  PRINT  S.N 

420.  5  FORMAT  (  UX,  *GH*  ,  12,  •  COOES  ARE:*) 
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421. 

PRINT  6.(01 1 > .1 =1 .MM) 

422. 

t 

FORMAT  (  '3  '  .  10* .  ‘HANDLES  ARE:  1 .61  ! 

423. 

70 

IF 1MM  . E0 .  N)  GO  TO  20 

424  . 

IF  (MM  .EO.  0>  GO  TO  10 

425. 

NUM  =  0 

426. 

IFL  =  0 

427. 

Ml  =  NPTS  *1 

426. 

1  *  1 

429. 

J  *  1 

430. 

DO  103  5  =  2. Ml 

431. 

IF  (5  .GT.  NPTS)  GO  TO  40 

432. 

FLAG  =  0 

433. 

11  =  ORDEROC) 

434. 

i 2  =  ORDER! K- 1 ) 

435. 

DO  200  Jl  =  1 .MM 

436. 

JJ  =  0(  Jl  ) 

437. 

!F(PTSET(  11  ,JJ>  .EO.  PTSETU, 

438. 

FLAG  =  l 

439. 

GO  TO  330 

440. 

200 

CONT ] NOE 

441  . 

330 

IF  (FLAG  .EO.  11  GO  TO  40 

442. 

J  *  J  ♦  1 

443. 

GO  TO  ICO 

444. 

40 

II  =  ORDER ( 1 1 

445. 

DO  300  Jl  *  l.MM 

446. 

JJ  =  0(JI  1 

447. 

C0DC1J1)  «  PT 5E 1(11. JJ) 

446. 

300 

CONT INOE 

449. 

Nl  =  MM* 1 

50. 

L  =  MM 

>1  . 

NN  =  N  -MM 

52. 

DO  400  Jl  *  N! . N 

33  • 

J  J  =  0(  Jl  ) 

454. 

CALL  M1NMAX1PTSET. order. JJ 

455. 

L  =  L  ♦  1 

456. 

CODE ( L  1  =  MIN 

457. 

400 

CDDEIL+NN)  k  MAX 

456. 

DEN  =  OCNST V ( CODE . MM . N. I .J) 

4S9. 

ifjoen  .ge. thresh)  go  to  m 

460. 

IFL  =  1 

461  . 

GO  TO  555 

462. 

1  1  1 

IF  ( I  .EO.  J 1  GO  TO  555 

463. 

N2  *  L 

4  64. 

N3  =  L  ♦  NN 

465. 

DO  333  15  =  I . J 

466. 

333 

ORDER ( 15  1=0 

467. 

IF  (  FLG  .EO.  1  )  GO  TO  72 

466. 

PRINT  32 , (CODE! 1 1 . I  «=  l.MM 

469. 

32 

TURMAT (  '0‘.  I  OX  .  1 (  ‘.613) 

470. 

PRINT  33 . (CODEC l 1 . 1*NI . N2 I 

471  . 

33 

FORMAT!  • ♦ • .30X. • ;  ‘.61  3) 

472. 

N22  =  N2»l 

4*3. 

PRINT  34, <CCDE< I 1 . I=N22 ,N3 

474  . 

34 

FORMAT (•♦•,50X.‘;‘.613> 

475. 

PRINT  35 

476. 

35 

FORMAT!  •  »60X» •)  • ) 

477. 

72 

NUM  *  NUM  «  | 

478. 

555 

I  =  5 

479. 

J  =  5 

480. 

1  OC 

CONT 1NUE 

GO  TO  200 


30 


rnmmm 


mm 
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*92. 

*93. 

*9*. 

*95. 

*95. 

*97. 

*99. 

•  99. 
490. 
*91  . 
*92  . 
*93. 
*9*. 

•  95. 
*95. 

•  97. 
*99. 
*99. 
500. 
501  . 

502. 

503. 
50*. 
505. 
505. 
507. 
509. 
**9. 

7. 


CCCCCCCCCC 

c 

c 

SUa°3UT!*e  EXHAUSIPTSFT.M.N.NCODE.IFCl 

I NTEGCR  PTSEM 110.61.0(6 I. TABLE ( 720.61 .NCODEI5. 31 

CALL  “RCnMBCTABLE.Nl 

ITER*  *  1 

no  loo  i  b  i .N 

100  l TERM  =  (TERM* I 

00  200  t  ■  I . I  TERM 
03  300  Jb  1 ,N 

300  0(J1  *  TABLE( 1 .Jl 

1F( IEL  .EO.  1  I  GO  TO  SB 
PRINT  l 

1  FORMAT (//. 1 IX , 'HANDLES  CHOSEN  ARE! » I 
PRINT  2.(0( Jl . J»1 .N > 

2  F3  4M4  T I 'O' « 10X.AI51 

99  CALL  SORT“T (OTSET .O.M.Nl 

|F  C IFL  .EO.  t  I  60  TO  89 
PRINT  3 

3  FORMAT!//, I  l  X. 'FOLLOWING  ARE  GH  COOES  W.R.T.  ABOVE  • 

•  .'HANDLES:'! 

89  N1  b  N  ♦  1 

03  500  J  =  1 . Nt 

C  ALL  ENCOOEi J.M.N.O.PTSET.NUM . IFL1 
I F ( J  .LE.I  .OR.  J  .GE.NII  GO  TO  500 
I F ( NUM  . L T .  NCODEI J-I.21I  NC00E1J-I.2I  »  NON 
TFIK'UN  .GT.  NCOOEI j-1,311  NCODEIJ-1,31  *  NUM 
500  CONTINUE 

!E(  IEL  .EO.  II  GO  TO  200 


-.1.  PRINT  A 

512.  *  FORMAT  I  /  /  ,  I  00 (*.'11 

513.  200  CONTINUE 

51*.  RETURN 

515.  END 

516.  CCCCCCCCCC 

517.  C 

518.  r  PRCOMBI T49LE.N1 :  t0  FIND  ALL  POSSIBLE  SET  OF  HANDLES.  USING 

519.  C  PERMUTATION  C OHBT N AT  I ON.E ACH  Row  OF  TABLE  CONTAINS  OK  SET 

520.  C  HANDLES. N  IS  THE  DIMENSION  OF  SPACE. 

521.  C  calls: 

522.  C  III  fill:  to  STORE  PROPER  COLUMN  numberes  in  table. 

5?3»  C 

524.  CCCCCCCCCC 

525.  SUBROUTINE  PRCOMBI T ABLE .N I 

526.  INTEGER  TABLET720.6I 

527.  INTEGER  COL 

528.  NUM  b  N 

529.  03  100  I  *  I.N 

530.  100  TABLE( t .1  I  b  I 

531.  K  b  | 

532.  COL  b  1 

533 •  10  L  =  I 

534.  N1  b  NUM  -  | 

535.  'TERM  b  1 

536.  f>0  203  I  =  |  ,NI 

537.  200  ITBRM  *  ITERM  »l 

538.  DO  300  Kl  *  t.K 

539.  03  409  I  b  I. NUM 

540.  I  I  b  L 
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